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ABSTRACT 

We have observed the OH (1720 MHz) line in five galactic SNRs with the VLA to 
measure their magnetic field strengths using the Zeeman effect. We detected all 12 of the 
bright (S u > 200 mJy) OH (1720 MHz) masers previously detected by Frail et al. (1996) 
and Green et al. (1997) and measured significant magnetic fields (i.e. > 3a) in ten of them. 
Assuming that the "thermal" Zeeman equation can be used to estimate | B | for OH masers, 
our estimated fields range from 0.2 to 2 mG. These magnetic field strengths are consistent 
with the hypothesis that ambient molecular cloud magnetic fields are compressed via the 
SNR shock to the observed values. Magnetic fields of this magnitude exert a considerable 
influence on the properties of the cloud with the magnetic pressures (10~ 7 — 10 -9 erg cm -3 ) 
exceeding the pressure in the ISM or even the thermal pressure of the hot gas interior to 
the remnant. This study brings the number of galactic SNRs with OH (1720 MHz) Zeeman 
detections to ten. 

Subject headings: 

ISMxlouds — ISM:individual (W51, G349.7+0.2, CTB37A, CTB33, G357.7-0.1) — ISM:magnetic 
fields — masers — polarization — radio lines:ISM 
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1. INTRODUCTION 

Magnetic fields can moderate the impact that 
a shock has on a molecular cloud. In the absence 
of a field, a supernova blast wave will heat, com- 
press and fragment the cloud and may ultimately 
destroy it (Klein, McKee & Colella 1994). The 
inclusion of a field mitigates these effects, lim- 
iting compression and stabilizing it against frag- 
mentation, allowing the cloud to survive, perhaps 
to trigger a future generation of star formation 
(Miesch & Zweibel 1994; MacLow et al. 1994). 

It has long been known that supernova rem- 
nants (SNRs) possess magnetic fields. Observa- 
tions of synchrotron radiation have established 
that the direction of the magnetic field in young 
SNRs like Cas A is predominately parallel to the 
shock normal (i.e. radial), whereas for older rem- 
nants the fields are perpendicular (e.g. Dickel et 
al. 1991; Milne 1990). The latter geometry likely 
originates from the compression of the ambient 
interstellar field, while Rayleigh- Taylor instabil- 
ities are invoked to explain the radial fields in 
young SNRs (Jun & Norman 1996). Until re- 
cently, estimates of the strength of the magnetic 
fields in SNRs had to rely on the somewhat du- 
bious equipartition approximation. This situa- 
tion has changed with the re-discovery of shock- 
excited OH (1720 MHz) maser emission in SNRs 
(Frail, Goss & Slysh 1994). 

In a series of recent papers, the satellite line of 
the OH molecule at 1720.53 MHz has been used 
as a powerful probe of SNR-molecular cloud in- 
teractions. OH (1720 MHz) masers are found 
in ~20 SNRs, or 10% of the known SNRs in 
our Galaxy (see Koralesky et al. 1998 and ref- 
erences therein). This maser line is inverted 
through collisions with H2 (n ~ fewxlO 4 cm~ 3 
and T~80 K) behind C-type shocks propagating 
into molecular clouds (Reach & Rho 1999, Frail 
& Mitchell 1998). The geometry of the shock is 
well constrained since strong maser amplification 
can only occur when the shock front is viewed 
edge on (Claussen et al. 1997). These obser- 
vational statements are well-supported by theo- 
retical modeling of the OH (1720 MHz) excita- 
tion (Lockett, Gauthier & Elitzur 1998; Wardle, 



Yusef-Zadeh, & Geballe 1999; Wardle 1999). 

One advantage of observing the OH (1720 
MHz) maser line is that it allows for a measure- 
ment of the strength of the magnetic field via 
Zeeman splitting (e.g. Troland & Heiles 1982). 
In this case, when the observed splitting is small 
compared to the line width, V = ZC\ B | dl/dz^, 
where | B | is the the total magnetic field 
strength, Z is the Zeeman splitting coefficient, 
and C is a constant which depends on the an- 
gle between the line of sight and B (the possible 
forms of C will be discussed in §4.1). For now we 
will denote the combination of C\ B | as Bg. 

This method has already been used to suc- 
cessfully measure the magnetic field strength in 
the post-shock gas behind five SNRs (Sgr A East, 
W44, W28, G32.8-0.1, and G346.6-0.2), yield- 
ing values for Bq between 0.1 to 4 mG (Yusef- 
Zadeh et al. 1996, Claussen et al. 1997, Ko- 
ralesky et al. 1998). With these promising re- 
sults in mind, we have performed 1720 MHz Zee- 
man studies toward five more of the 17 SNRs 
found in surveys by Frail et al. (1996) and Green 
et al. (1997) to contain OH (1720 MHz) masers 
using the NRACf] VLA. Only SNRs with bright 
masers (I>200 mJy) and no previous high reso- 
lution Stokes V observations were chosen. 

2. OBSERVATIONS 

We observed five SNRs (W51, G349.7+0.2, 
CTB37A, CTB33, and G357.7-0.1) at 1720 MHz 
with the VLA in A configuration. The key ob- 
serving parameters for each SNR are given in 
Table 1. All of the SNRs were observed in 
"2IF mode" (recording both right (RCP) and left 
(LCP) circular polarization) with a 0.1953 MHz 
(34.0 km s _1 ) bandwidth divided into 127 chan- 
nels. Due to this narrow bandwidth, CTB37A 
had to be observed with two different center fre- 
quencies (-22 km s _1 and -64 km s _1 ) in order to 
obtain data for its two different OH (1720 MHz) 
maser populations. The observations were con- 

3 The National Radio Astronomy Observatory is a facil- 
ity of the National Science Foundation operated under 
a cooperative agreement by Associated Universities, 
Inc. 
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Table 1 

Observational Parameters with the VLA 



SNR 


R. A. a 


Decl. a 


Vi sr 


t source 


Beam 


P. A. 


a b 

u rms 




(B1950) 


(B1950) 


(km s _1 ) 


(hours) 


(" x ") 


(°) 


(mjy) 


W51C (2IF) 


19 20 35.0 


+ 14 09 00.0 


+64.0 


3.3 


1.4 x 1.2 


-40 


6.6 


(PA) C 


19 20 35.0 


+ 14 09 00.0 


+70.0 


1.6 


1.6 x 1.3 


+55 


6.5 


G349.7+0.2 


17 14 36.0 


-37 22 00.0 


+16.0 


2.6 


3.4 x 1.3 


+08 


7.2 


CTB37A (I) 


17 11 00.0 


-38 30 00.0 


-64.0 


4.9 


3.9 x 1.2 


-14 


5.4 


(II) 


17 11 00.0 


-38 30 00.0 


-22.0 


4.0 


4.2 x 2.3 


+23 


5.5 


CTB33 


16 32 14.0 


-47 31 00.0 


-70.0 


2.1 


8.0 x 1.1 


+67 


10.2 


G357.7-0.1 


17 36 56.0 


-30 56 00.0 


-12.0 


5.1 


3.0 x 1.2 


+ 18 


5.2 



a Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. 

b rms noise in an individual channel. 

C W51C (PA) data have a velocity resolution of 0.56 km s _1 while all other data have velocity resolution of 
0.27 km s" 1 . 



ducted on June 22, July 6, 8, 12, 15, September 
13, and October 24, 1999. The OH data were 
Hanning smoothed online, and the resulting ve- 
locity resolution is 0.27 km s -1 . We observed 
both senses of circular polarization simultane- 
ously and since Zeeman observations are very 
sensitive to small variations in the bandpass, a 
front-end transfer switch was used to periodically 
switch the sense of circular polarization passing 
through each telescope's IF system. 

The AIPS (Astronomical Image Processing 
System) package of the NRAO was used for the 
calibration, imaging, and cleaning of the OH 
(1720 MHz) data sets. The RCP and LCP 
data were calibrated separately and later com- 
bined during the imaging process to make Stokes 
I=(RCP + LCP)/2 and Stokes V=(RCP - LCP)/2 
data sets. Bandpass correction was applied only 
to the I data sets since bandpass effects cancel 
to first order in the V data. Line data sets were 
created by estimating and removing the contin- 
uum emission in the UV plane using the AIPS 
task UVLSF. The strongest maser channel in 
each line data set was then self-calibrated and 
the solutions were applied to each channel. Sub- 



sequent magnetic field estimates were performed 
using the MIRIAD (Multichannel Image Recon- 
struction Image Analysis and Display) processing 
package from BIMA. The rms noise per spectral 
channel obtained for each SNR is summarized in 
Table 1. 

After imaging, each of the bright maser spots 
was fit with a 2-D Gaussian using the AIPS task 
JMFIT. None of these fits showed convincing ev- 
idence that the individual maser spots are re- 
solved at the resolutions shown in Table 1; there- 
fore, we regard the maximum sizes reported by 
JMFIT to be upper limits. The positions, peak 
flux densities, and an upper limit to the maser 
spot sizes are reported in Table 2. In addi- 
tion, the Stokes I spectrum at the peak pixel of 
each maser spot was fit with a single Gaussian 
in the spectral domain using GIPSY, to obtain 
each maser's center velocity (t>) and linewidth 
(Avfwhm)- These estimates of Avfwhm were 
then corrected for the finite channel width of 
the the data (0.27 km s _1 ). These deconvolved 
Avfwhm and center velocities are also reported 
in Table 2. The absolute position errors of these 
data are ~ 0.1", while the relative position er- 
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Table 2 

Fitted Parameters of OH (1720 MHz) Masers 



SNR 


Feature 


R. A. a 


Decl. a 


C b 

>~>peak 


Vlsr ° 


Avfwhm d 


$max 






(B1950) 


(B1950) 


(mjy) 


(km s _1 ) 


(km s _1 ) 


(") 


W51C 


1 


19 20 35.9 


+14 09 53.6 


2710 


+72.0 


0.9 


0.7 




2 


19 20 36.4 


+14 09 50.3 


4760 


+69.1 


1.2 


0.1 


G349.7+0.2 .. 


1 


17 14 36.0 


-37 23 01.9 


770 


+16.2 


1.1 


1.6 




2 


17 14 36.9 


-37 22 52.8 


1800 


+15.2 


0.3 


1.0 




3 


17 14 37.5 


-37 23 14.7 


1060 


+16.7 


0.6 


0.9 


CTB37A 


1 


17 10 51.8 


-38 28 51.7 


410 


-66.2 


0.7 


0.6 




2 


17 10 56.8 


-38 28 56.7 


730 


-63.7 


0.5 


0.3 




3 


17 10 59.3 


-38 31 19.2 


470 


-63.5 


1.5 


0.2 




4 


17 11 01.4 


-38 37 32.9 


220 


-65.3 


1.1 


0.7 




5 


17 11 09.0 


-38 25 51.3 


400 


-21.5 


0.9 


1.7 


CTB33 


1 


16 32 06.2 


-47 29 53.2 


250 


-71.8 


1.0 


1.7 


G357.7-0.1 .. 


1 


17 36 54.6 


-30 56 07.6 


400 


-12.3 


0.6 


0.8 



a Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. 

b Errors in the peak flux range from 4 to 8 mjy. 
c Errors in vi sr range from 0.02 to 0.002 km s _1 . 

d Errors in Avfwhm range from 0.05 to 0.003 km s _1 , and the values shown have been deconvolved from 
the finite channel width (0.27 km s _1 ). 
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rors (compared to other masers in the field) are 
~ 0.02". 

In addition, W51 was observed on August 3, 
1999 in "PA mode" (recording all Stokes parame- 
ters) with a 0.1953 MHz bandwidth, and 64 chan- 
nels. This correlator setup resulted in a veloc- 
ity resolution of ~ 0.54 km s _1 after Hanning 
smoothing. The details of this observation can 
also be found in Table 1. These data were cali- 
brated in the same manner described above, with 
the exception of the polarization calibration. The 
absolute polarization calibration was carried out 
by extrapolating 3C48 data from seven previous 
20 cm polarimetry observations to yield a posi- 
tion angle of P.A.= 13° for this calibrator at 1720 
MHz. The error in this estimate (~ 5°) will also 
apply to the position angles derived from these 
data. 

3. RESULTS 

3.1. General Maser Properties 

From the five observed SNRs, we detected 
all 12 bright OH (1720 MHz) masers (S > 200 
mJy) previously known from the surveys of Frail 
et al. (1996) and Green et al. (1997). Ta- 
ble 2 summarizes the fitted parameters of these 
maser features. Of these, ten show classical S- 
shaped Stokes V profiles (e.g. Elitzur 1998), 
one shows a complicated Stokes V profile, and 
one shows no discernible V signal despite its 
~ 1800 mJy peak flux density. The observed 
masers have typical deconvolved line widths of 
~ 0.8 km s _1 , spanning the range from 0.3 to 
1.5 km s _1 . In addition, none of the observed 
OH (1720 MHz) masers appear to have under- 
gone significant changes in flux density or po- 
sition since the detection experiments were per- 
formed in 1994. As noted in §2, none of the maser 
spots have been resolved at the resolutions listed 
in Table 1. This may be due, in part, to the 
low declination, and hence elliptical synthesized 
beam shapes observed for four of the five SNRs. 
Using the upper limits on the maser spot sizes 
listed in Table 2, lower limits on the brightness 
temperatures of these masers range from 4 x 10 4 
K to 10 8 K. 



Using the thermal Zeeman expression (V = 
ZC\ B | dl/du) with Z = 0.6536 Hz ^G" 1 for 
OH at 1720 MHz, the derivative of Stokes I was 
fitted to Stokes V in a least square fitting rou- 
tine as described by Claussen et al (1997), to 
obtain the combination C\ B \=Bq. The validity 
of the thermal Zeeman formulation for the mag- 
netic field in masers along with the 9 dependence 
and magnitude of C is discussed in §4.1. Mag- 
netic fields were detected toward all five SNR at 
greater than the 3a level. The fitted magnetic 
field strengths (Be), and their associated la er- 
rors are summarized in Table 3, while comments 
on individual sources appear in §3.2. The de- 
tected Bg fields range from 0.2 to 2 mG. Un- 
like W28 and W44, where the fields were found 
to be uniform in both magnitude and direction 
(Claussen et al. 1997), CTB37A shows a compli- 
cated magnetic field morphology. In this source 
Bg changes by a factor of seven and reverses sign 
on a length scale of ~ 3 pc. 

3.2. Individual Sources 

3.2.1. W51 

W51 is composed of two H n region complexes 
(W51A and W51B), as well as the SNR W51C 
and is located at the tangent point of the Sagit- 
tarius arm at I = 49° corresponding to a dis- 
tance of ~ 6 kpc (Koo, Kim, & Seward 1995). A 
continuum image of W51 at 330 MHz from Sub- 
rahmanyan & Goss (1995) is shown in Figure 1 
with the two W51 OH (1720 MHz) masers super- 
posed. The spectral index of W51C is difficult to 
calculate due to contamination from the W51B 
Hll regions, but Subrahmanyan & Goss (1995) 
estimate that a = —0.2 (S = v a ). 

Based on the absorption of X-ray emission 
from W51C toward W51B, Koo et al. (1995) sug- 
gest that W51C must be located behind W51B. 
Koo k Moon (1997a, 1997b) observed high ve- 
locity (HV) Hi, CO(1-0), and CO(2-l) emis- 
sion between +85 and +120 km s _1 to the east 
of W51B where it overlaps W51C (see Fig. 1). 
They interpret this HV gas as arising from a 
shock interaction between W51C and a molecular 
cloud (located between W51C and W51B, or pos- 



5 



Table 3 

Magnetic Fields 



SNR 


Feature 


B e a 






(mG) 


W51C 


1 


+ 1.5 ± 0.05 




2 


+ 1.9 ± 0.10 


G349.7+0.2 ... 


1 


complex 6 




2 


<;o.i c 




3 


-0.35 ± 0.05 


CTB37A 


1 


-0.5 ± 0.10 




2 


+0.22 ± 0.05 




3 


-0.60 ± 0.09 




4 


+ 1.5 ± 0.20 




5 


-0.8 ± 0.10 


CTB33 


1 


+ 1.1 ± 0.30 


G357.7-0.1 .... 


1 


+0.7 ± 0.12 



a These Be field estimates were calculated using the 
thermal Zeeman equation and may overestimate the 
maser fields according to the maser theory of Elitzur 
(1998). The correction cannot be determined exactly 
but is probably on the order of (0.5 — 0.2)xBg, there- 
fore these values should be viewed as upper limits. See 
§4.1 for further discussion. 

b Stokes V profile is complex, probably indicative of 
blending. 

c Upper limit to Bg based on S/N of data. See also 
§3.2.2. 
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19 h 22 m 30 s 21 m 30 s 20 m 30 s 1 9 m 30 s 
Right Ascension (1950) 



Fig. 1. — W51 VLA 330 MHz continuum image 
from Subrahmanyan & Goss (1995) with contours at 
(0.02, 0.12, 0.22, 0.32, 0.42, 0.52, 0.62, 0.72, 0.82) 
x 2.46 Jy beam -1 . The resolution of this image is 
- V. The locations of the two W51C OH (1720 MHz) 
masers are indicated by the white (+) symbols. 

sibly the backside of the W51B complex itself). 
They show that the HV H I is located toward the 
western edge of the centrally bright X-ray emit- 
ting region of W51C in an arc-shape, while the 
shocked CO gas is located slightly east of the HV 
Hi emission. The two OH (1720 MHz) masers 
reported in this work (and Green et al. 1997) 
are located toward the NE end of the shocked 
CO/Hi arc structure and ~ 2' SW of the Hn re- 
gion G49.2-0.3 (see Fig. 1) at velocities of ~ 70 
km s _1 . From the coincidence of these OH (1720 
MHz) masers with the shocked gas and their an- 
gular separation from G49.2— 0.3 (~ 3.5 pc for 
d = 6 kpc), Green et al. (1997) argue that they 
are associated with the W51C SNR shock. 

Fits of B for the two W51C OH (1720 MHz) 
masers are shown in Figure 2. The values of Bg 
for these masers are 1.5 and 1.9 mG for features 
1 and 2, respectively (see Tables 2 and 3). The 
W51C masers were also observed in PA mode, 
providing images of all four Stokes parameters (I, 
V, Q, U), with 0.54 km s -1 velocity resolution. 
Toward the strongest W51C OH maser (feature 
2), Stokes Q and U signals were detected at 15<r 




»UR ft™ ■"') "ls. 0»n •"') 



Fig. 2. — Fits of B e for W51 OH (1720 MHz) maser 
features (1) and (2). The upper panels show the VLA 
Stokes I profiles (solid histogram), and the bottom 
panels show the VLA Stokes V profiles (solid his- 
togram) with the fitted derivative of Stokes I shown 
as smooth dotted curves. The value of B e fit for each 
position and its calculated error are given at the top 
of each plot. 

and 17a, respectively. A profile showing Stokes 
Q and U toward this maser is shown in Figure 
3. Using Stokes Q and U values of Q = 0.10 Jy 
beam -1 and U = —0.11 Jy beam -1 , we obtain a 
linear polarization percentage of 3.5% and a P.A. 
of —25° for the linear polarization vector. 

3.2.2. G349.7+0.2 

G349.7+0.2 is one of the most luminous SNR's 
in the galaxy (after Cas A and the Crab; Shaver 
et al. 1985a), if it is located at a distance of 
~ 22 kpc (Frail et al. 1996). The spectral in- 
dex of G349.7+0.2 was estimated by Shaver et al. 
(1985a) to be ~ —0.5, typical of shell type rem- 
nants. This SNR contains three bright maser fea- 
tures along with several weaker features within 
~ 1' (6 pc) near the center of the remnant. The 
positions of the brightest masers are shown on 
a continuum image of G349.7+0.2 in Figure 4. 
This continuum image is a compilation of 18 and 
20 cm data from the VLA archive, and contains 
data from all four configurations (A, B, C, and 
D). It is the most sensitive and highest resolu- 
tion continuum image of this SNR to date, with 
a resolution of 5" x 2" (P.A.= 0.6°) and rms noise 
of ~ 1 mJy beam -1 . 
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Fig. 3. — VLA Stokes Q and U linear polarization 
profiles (upper panel) and Stokes V circular polariza- 
tion profile {lower panel) for W51C OH (f 720 MHz) 
maser feature (2). Note that the spectral resolution 
of these data are only 0.56 km s _1 while the spec- 
tra displayed in Fig. 2 (and all other Figures) have a 
resolution of 0.27 km s -1 . 

Only G349.7+0.2 OH(3) (OH maser feature 
3) has a significant magnetic field detection with 
Bq= 0.35 mG. The Bg fit for this feature is shown 
in Figure 5. Maser features 1 and 2 are the only 
bright (S > 200 mJy) OH (1720 MHz) masers in 
our sample for which we were unable to detect 
a significant B e . G349.7+0.2 OH(l) exhibits a 
complicated Stokes V spectrum that is indicative 
of blending. Therefore, it is possible that im- 
proved velocity and/or angular resolution would 
lead to a Bg detection for this maser. G349.7+0.2 
OH(2) shows no hint of any Stokes V signal de- 
spite its 1800 mJy flux density. The S/N of these 
data alone suggest that Bq< 0.1 mG, however, 
this maser is also the narrowest maser in our 
sample (Avfwhm = 0.3 km s -1 ). Therefore, it 
is possible that if the Stokes V was not resolved 
in velocity, B e for G349.7+0.2 OH(2) could be 
~ 0.3 mG. 

3.2.3. CTB37A 

CTB37A, also known as G348.5+0.1, is esti- 
mated to lie at a distance of ~ 11 kpc (Frail et al. 
1996). Kassim, Baum, k, Weiler (1991) estimate 
its spectral index to be ~ —0.5 based on flux 
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Fig. 4. — G349.7+0.2 20 cm VLA continuum image 
with contours at (0.05, 0.15, 0.25, 0.35, 0.45, 0.55, 
0.65, 0.75, 0.85, 0.95) x 63 mJy beam" 1 . This image 
is composed of archival VLA data from A, B, C, and 
D configurations. The resolution is 5" x 2" (P.A. = 
0.6°) and the peak flux density is 63 mJy beam -1 
with an rms noise of ~ 1 mJy beam -1 . The locations 
of the three G349.7+0.2 OH (1720 MHz) masers with 
S v > 200 mJy are indicated by the black (+) symbols. 

density measurements ranging from 80 MHz to 
14.7 GHz. A 21 cm continuum map of CTB37A 
from Kassim et al. (1991) is displayed in Fig- 
ure 6 with the OH (1720 MHz) masers super- 
posed. CTB37A contains two kinematically dis- 
tinct sets of OH (1720 MHz) masers. One group 
at ~ —22 km s -1 is located toward the north 
end of CTB37A, while the others have velocities 
of ~ —65 km s _1 and are located near the center 
and southern parts of the source. Kassim et al. 
(1991) propose that the extension of continuum 
emission seen to the east of CTB37A (see Fig. 6) 
is a separate SNR which they name G348.5— 0.0. 
It was further proposed by Frail et al. (1996) 
that the ~ —22 km s -1 masers originate from 
this second SNR. 

CTB37A, also known as G348.5+0.1, is esti- 
mated to lie at a distance of ~ 11 kpc (Frail et al. 
1996). Kassim, Baum, &: Weiler (1991) estimate 
its spectral index to be ~ —0.5 based on flux 
density measurements ranging from 80 MHz to 
14.7 GHz. A 21 cm continuum map of CTB37A 
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Fig. 5.— Fit of B e for G349. 7+0.2 OH (1720 MHz) 
maser feature (3). The upper panels show the VLA 
Stokes I profiles [solid histogram), and the bottom 
panels show the VLA Stokes V profiles (solid his- 
togram) with the fitted derivative of Stokes I shown 
as smooth dotted curves. The value of Be fit for each 
position and its calculated error are given at the top 
of each plot. 

from Kassim et al. (1991) is displayed in Fig- 
ure 6 with the OH (1720 MHz) masers super- 
posed. CTB37A contains two kinematically dis- 
tinct sets of OH (1720 MHz) masers. One group 
at ~ —22 km s^ 1 is located toward the north 
end of CTB37A, while the others have velocities 
of ~ —65 km s _1 and are located near the center 
and southern parts of the source. Kassim et al. 
(1991) propose that the extension of continuum 
emission seen to the east of CTB37A (see Fig. 6) 
is a separate SNR which they name G348.5— 0.0. 
It was further proposed by Frail et al. (1996) 
that the ~ —22 km s _1 masers originate from 
this second SNR. 

As noted at the beginning of this section, the 

B e of the 65 km s" 1 OH (1720 MHz) masers 

have the unusual property that they change di- 
rection along the line of site over length scales 
as small as 1' (~ 3 pc). In addition, the magni- 
tude of Bq changes by a factor of seven between 
CTB37A OH(2) and CTB37A OH(4) from ~ 0.2 
mG to ~ 1.5 mG. Fits of Bq for these four ~ —65 
km s _1 masers are shown in Figure 7. The only 
strong ~ —22 km s -1 OH maser has Bq= —0.8 
mG and the fit for this maser feature (CTB37A 
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Fig. 6. — CTB37A 20 cm VLA continuum image 
from Kassim et al. (1991) with contours at (0.01, 
0.12, 0.23, 0.34, 0.45, 0.56, 0.66, 0.76 0.87) x 0.427 
Jy bcaur 1 . The beam is 33" x 18" (P.A.= 5°). The 
(+) symbols mark the locations of the ~ —65 km s _1 
OH (1720 MHz) maser features, while the locations 

of the 22 km s" 1 OH (1720 MHz) maser features 

are marked with (x) symbols. 

OH(5)) is also shown in Fig. 7. 

Reynoso & Mangum (1999, in prep.) have de- 
tected CO(1-0) emission toward CTB37A with 
~ 1' resolution using the Kitt Peak 12m tele- 
scope. Their maps show two distinct CO clouds, 
one at ~ —22 km s _1 in the northern part of 
CTB37A, and another at ~ —65 km s _1 which is 
concentrated to the northwest and middle of the 
source. Both CO clouds are coincident spatially 
and in velocity with our two groups of OH (1720 
MHz) masers. 

As noted at the beginning of this section, the 

B e of the 65 km s" 1 OH (1720 MHz) masers 

have the unusual property that they change di- 
rection along the line of site over length scales 
as small as 1' (~ 3 pc). In addition, the magni- 
tude of Bq changes by a factor of seven between 
CTB37A OH(2) and CTB37A OH(4) from ~ 0.2 
mG to ~ 1.5 mG. Fits of Bq for these four ~ —65 
km s _1 masers are shown in Figure 7. The only 
strong ~ —22 km s _1 OH maser has Bq= —0.8 
mG and the fit for this maser feature (CTB37A 
OH(5)) is also shown in Fig. 7. 
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Fig. 7.— Fits of B e for CTB37A OH (1720 MHz) 
maser features (1 - 5). Note that features (1-4) cor- 
respond to the ~ —65 km s^ 1 masers, while feature 
(5) is the single strong <~ —22 km s _1 maser. The up- 
per panels show the VLA Stokes I profiles (solid his- 
togram), and the bottom panels show the VLA Stokes 
V profiles (solid histogram) with the fitted derivative 
of Stokes I shown as smooth dotted curves. The value 
of Bq fit for each position and its calculated error are 
given at the top of each plot. 



Reynoso &; Mangum (1999, in prep.) have de- 
tected CO(l-O) emission toward CTB37A with 
~ 1' resolution using the Kitt Peak 12m tele- 
scope. Their maps show two distinct CO clouds, 
one at ~ —22 km s _1 in the northern part of 
CTB37A, and another at ~ —65 km s™ 1 which is 
concentrated to the northwest and middle of the 
source. Both CO clouds are coincident spatially 
and in velocity with our two groups of OH (1720 
MHz) masers. 

3.2.4. CTB33 

The SNR in CTB33 is also known as G337.0-0.1, 
and was estimated to lie at a distance of ~ 11 




Fig. 7. — Continued... 

kpc by Sarma et al. (1997) using Hi absorption 
data from the ATCA (Australia Telescope Com- 
pact Array) . These authors also find that CTB33 
has a spectral index of ~ —0.6, typical of shell 
type SNRs. The locations of the two bright OH 
(1720 MHz) masers from Frail et al. (1996) are 
shown on Figure 8 superposed on a 1380 MHz 
CTB33 continuum image (Sarma et al. 1997). 
However, Sarma et al. (1997) present evidence 
that the southern most CTB33 maser is proba- 
bly associated with the Hn region G336.9— 0.2. 
Our OH (1720 MHz) results are consistent with 
this suggestion since this maser feature has been 
resolved into at least four different velocity com- 
ponents, each of which show complicated Stokes 
V spectra typical of Hn region masers (i.e. no 
"S" shaped Zeeman pattern; see Elitzur 1998). 
For this reason only the centrally located CTB33 
maser (CTB33 OH(l)) has been included in Ta- 
bles 2 and 3. The Bg detected for this maser is 1.1 
mG, but the range of channels over which the fit 
was performed had to be limited to obtain a rea- 
sonable fit (see Fig. 9). Therefore, this detection 
should be considered tentative. Such a suppres- 
sion of part of the Stokes V signal may be the re- 
sult of blending of the maser line in velocity space 
since this maser is fairly wide (Avfwhm^ 1 
km s^ 1 , also see Roberts et al. 1993), but could 
also simply be a result of poor S/N. 

Corbel et al. (1999) have detected a number of 
molecular clouds toward CTB33 in CO( J = 1—0) 
and CO(J = 2-1) emission using the SEST tele- 
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Right Ascension (1950) 

Fig. 8. — CTB33 1375 MHz ATCA continuum im- 
age from Sarma et al. (1997) with contours at (0.3, 
0.4, 0.5, 0.6, 0.7, 0.8, 0.9) x 0.090 Jy beam" 1 . The 
resolution of this image is <~ 12". The locations of the 
two CTB33 OH (1720 MHz) masers with S v > 200 
mJy are indicated by the black (+) symbols. 

scope. In particular, they find a giant molecular 
cloud at a velocity of ~ —71 km s _1 (Av = 11 
km s _1 ) with an approximate size of 67 pc and a 
mass of 4 x 10 6 Mq. The velocity of this molec- 
ular cloud is in excellent agreement with the ve- 
locity of CTB33 OH(l) (-71.8 km s' 1 ) and may 
well be the origin of the SNR/molecular cloud 
interaction needed to pump the maser (see §1 & 
§4.2). These authors also suggest that the SNR 
must lie on the near-side of the ~ —71 km s _1 
GMC based on comparison of their CO data (and 
its implied extinction) with the X-ray data from 
Woods et al. (1999). 

The CTB33 SNR is coincident with the site of 
a recently discovered Soft Gamma ray Repeater 
(SGR1627-41; see Hurley et al. 1999). The loca- 
tion of the SGR has been constrained by the 3rd 
Interplanetary Network (IPN: Ulysses, KONUS- 
WIND and BATSE) along with the detection of 
a presumably related BeppoSAX X-ray source to 
lie within the error boxes shown in Figure 10 on 
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Fig. 9.— Fit of B e for CTB33 OH (1720 MHz) 
maser feature (1). The upper panels show the VLA 
Stokes I profiles {solid histogram), and the bottom 
panels show the VLA Stokes V profiles {solid his- 
togram) with the fitted derivative of Stokes I shown 
as smooth dotted curves. The value of Bq fit for each 
position and its calculated error are given at the top 
of each plot. The solid portion of the Stokes I his- 
togram {upper panel) shows the velocity range used 
in the fit. 

a simplified contour image of CTB33 (Sarma et 
al. 1997; Hurley et al. 1999; Woods et al. 1999). 
Current theories for the nature of SGR's suggest 
that they arise from strongly magnetized neutron 
stars or 'magnetars', and that the outbursts are 
the result of crustquakes on the surface of the 
neutron star (Thompson & Duncan 1995). The 
discovery of SGR1627— 41, marks only the fourth 
such source to be detected. 

All of the previously known SGRs have been 
associated with young SNRs (although SGR1900+14 
lies close to, but not inside its associated SNR; 
see Hurley et al. 1999 and references therein). 
For this reason, despite the inexact correspon- 
dence between the IPN/BeppoSAX error boxes 
and the extent of CTB33 (95") proposed by 
Sarma et al. (1997), SGR1627-41 has been 
assumed to be the progenitor of CTB33. For 
this reason, a number of authors (Hurley et al. 
1999; Corbel et al. 1999) have estimated that 
SGR1627— 41 must have a transverse velocity of 
more than ~ 1,000 km s _1 if CTB33's age is 

5, 000 years. Although these velocities are not 
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Right Ascension (J2000) 

Fig. 10. — Similar to Figure 8, but in J2000 co- 
ordinates with contours at (0.4, 0.6, 0.8) x 0.090 Jy 
beam -1 . The IPN and BeppoSAX error boxes are 
superposed showing the location of SGR1627— 41. 

unreasonable (Lyne & Lorimer 1994), the loca- 
tion of the maser CTB33 OH(l) suggests that 
CTB33 actually extends farther west as shown 
in Figure 10, and has a 'blowout' morphology 
toward the NE and SW. Such a morphology 
would not be unexpected in a region populated 
by H II regions that could have effectively cleared 
such cavities (see Jones et al. 1998). There- 
fore, an increase in the estimated size of CTB33 
would place the IPN/BeppoSAX error boxes and 
the OH (1720 MHz) maser much closer to the 
center of CTB33, rather than at its outskirts, 
with a substantial reduction in the estimate for 
the transverse velocity of SGR1627— 41. Future 
molecular, or high resolution, low-frequency con- 
tinuum data toward CTB33 may help resolve this 
issue. 

3.2.5. G357.7-0.1 (Tornado) 

G357.7-0.1 is an unusual SNR candidate lo- 
cated near the galactic center with a non-thermal 
spectral index in the range —0.5 < a < — 1.0 



(Stewart et al. 1994). It has been variously con- 
sidered to be everything from an extragalactic 
head-tail or double lobed source (Weiler & Pana- 
gia 1980; Caswell et al. 1989) to a new class of 
galactic head-tail object (Becker & Helfand 1985; 
Helfand & Becker 1985; see also Gray 1994). An- 
other intriguing suggestion is that G357.7— 0.1 is 
powered by an object ejected from the nearby 
SNR (G359.0-0.9) which lies only 1' from the 
symmetry axis of G357.7— 0.1. This scenario, 
however, would require such a 'runaway' pul- 
sar or X-ray binary to have a transverse veloc- 
ity ^ 2, 000 km s -1 which is considered unlikely 
(Gray 1994). The discovery of a OH (1720 MHz) 
maser coincident with G357.7— 0.1 has renewed 
speculation that it is, in fact, a galactic SNR 
(Frail et al. 1996). 

The odd morphology of G357.7— 0.1 has earned 
it the name 'Tornado'. The reason for this 
moniker can be seen in the greyscale continuum 
images of G357.7— 0.1 shown in Figure 11 (see 
also, Shaver et al. 1985a; Becker & Helfand 
1985). Figure 11a was created from archival VLA 
data at ~ 6 cm from both C and D configu- 
ration data. The resulting image has a resolu- 
tion of 13.4" x 6.7" (P. A. = 23.8°) and an rms 
of ~ 0.6 mJy beam -1 . The data used to create 
the image in Figure lib were compiled from 18 
cm and 20 cm VLA archive data from all four 
VLA configurations. The resulting image is one 
of the most sensitive and highest resolution im- 
ages of this source to date with an rms of ~ 1.5 
mJy beam -1 and resolution of 11". 6 x 8".l (P.A. 
= 26.2°; see also Shaver et al. 1985a; Becker 
& Helfand 1985; Yusef-Zadeh 1999b). Com- 
parison of the peak flux densities in Figs. 11a 
(£6cm =77 mJy beam -1 ) and lib (^ocm =145 
mJy beam -1 ) confirms that the spectral index of 
G357.7-0.1 is « -0.5. 
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Fig. 11. — Greyscale continuum images of 
G357.7— 0.1 composed of archival VLA data showing 
its spiral morphology, a) 6 cm image with a resolution 
of 13" A x 6". 7 (P.A. = 23.8°) and a peak flux density 
of 77 mJy beam -1 with an rms noise of ~ 0.6 mJy 
beam -1 , b) Same as a) but at 20 cm. The resolution 
of the 20 cm image is 11".6 x 8".l (P.A. = 26.2°) and 
the peak flux is 145 mJy beam -1 with an rms noise 
of ~ 1.5 mJy beam -1 . 
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Fig. 12. — Similar to Figure lib, but with contours 
at (0.05, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85) 
x 145 mJy beam -1 . The location of G357.7-0.1 OH 
(1720 MHz) maser feature (1) is marked by an (x) 
symbol. 



G357. 7-0.1 (1) 
= 690 ± 120 
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Fig. 13. — Fit of B e for G357.7-0.1 OH (1720 MHz) 
maser feature (1). The upper panels show the VLA 
Stokes I profiles (solid histogram), and the bottom 
panels show the VLA Stokes V profiles (solid his- 
togram) with the fitted derivative of Stokes I shown 
as smooth dotted curves. The value of Be fit for each 
position and its calculated error are given at the top 
of each plot. The solid portion of the Stokes I his- 
togram (upper panel) shows the velocity range used 
in the fit. 
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A compact source at the western edge of 
G357.7— 0.1 is also apparent in Figs. 11a and 
b (also observed by Shaver et al. (1985b); see 
also Stewart et al. 1994, and references therein). 
We obtain flux densities of 70 ±4 mJy and 69 ±12 
mJy for this compact source at 6 cm and 20 
cm respectively, confirming that this source has 
a flat spectrum as suggested by Shaver et al. 
1985b. At the resolution of the 18 cm A-array 
data (4". 5 x 3". 4) this source is extended and 
has a diameter of ~ 6" with a possible exten- 
sion to the east (toward the SNR). The exact 
nature of this source remains unknown. Shaver 
et al. (1985b) cite its flat spectral index at radio 
wavelengths as evidence that it is an unrelated 
Hn region. Alternatively, Shull, Fesen, & Saken 
(1989) suggest that it could be a pulsar which 
is interacting with the SNR shell and is respon- 
sible for the peculiar morphology of G357.7— 0.1 
similar to CTB 80. Higher resolution continuum 
observations or recombination line observations 
of this compact source are needed to determine 
its nature. 

Figure 12 shows the 20 cm image (similar 
to Fig. lib) with the single G357.7-0.1 OH 
(1720 MHz) maser superposed. The fit of this 
maser's Bg = 0.7 mG is shown in Figure 13. 
For this maser feature, only the high velocity 
side of the line could be adequately fit by the 
derivative of Stokes I. As mentioned for CTB33 
OH(l), a suppression of part of the Stokes V 
signal can be due to the presence of more than 
one velocity component, but the narrowness of 
this line (Avfwhm^ 0.6 km s _1 ) makes either 
an inability to fully resolve the Stokes V pat- 
tern or poor S/N more probable explanations. 
In the absence of better data, this detection 
should be regarded as tentative. Yusef-Zadeh et 
al. (1999b) also observed G357.7-0.1 at 1720 
MHz with the VLA in A-configuration during 
an OH (1720 MHz) survey of the Galactic cen- 
ter. With their short integration time (rms ~ 12 
mJy), they were able to place an upper limit 
on Bg for G357.7-0.1 OH(l) of 2 mG, consis- 
tent with our 0.7 mG detection. In addition, 
these authors detected extended OH (1720 MHz) 
emission toward G357.7— 0.1 with the VLA in D- 



configuration (114" x 38" resolution). These au- 
thors note that this extended emission may orig- 
inate from low-gain masers given the lack of OH 
(1720 MHz) absorption toward this source. 

Linear polarization images of the synchrotron 
emission from G357.7— 0.1 at 5.8 GHz by Stewart 
et al. (1994) show that the plane-of-sky magnetic 
field vectors lie circumferentially to the vertical 
bands of continuum emission (best seen in Fig. 
11). These authors note that this morphology is 
suggestive of a spiral magnetic field structure in 
the SNR. It is also interesting to note that the 
Stewart et al. images indicate a lack of linear po- 
larization at the location of the OH (1720 MHz) 
maser, as might be expected if B were nearly 
along the line of sight. 

4. DISCUSSION 

4.1. Implications from 1720 MHz Maser 
Theory 

We have assumed until now that magnetic 
fields can be calculated from the Zeeman effect in 
masers using the same formalism that has been 
developed for thermal radiation (i.e. Troland 
& Heiles 1982 and Bg in Table 3). However, 
stimulated emission is intrinsically different from 
spontaneous emission or absorption (see Elitzur 
1996), so it is not unreasonable to imagine that a 
maser's polarization properties are different from 
those of thermal emission. A common variable 
in any polarization study is the ratio of the line 
splitting (Aub = Z \ B |) to the Doppler width 
(Aud) which we will denote xb- In what fol- 
lows we assume this ratio is small {xb « 1) 
since we do not resolve the line splitting in our 
Stokes I spectra (i.e. the splitting is only de- 
tected in Stokes V) . In this section we review the 
results of several theoretical maser polarization 
studies for the case xb « 1 (Elitzur 1996, 1998; 
and Nedoluha & Watson 1992) and compare the 
implications and predictions of these studies to 
the currently available OH (1720 MHz) Zeeman 
data. These comparisons have not proven suffi- 
cient to discern exactly how the total maser mag- 
netic field strength | B | scales with our thermal 
estimates Bg (Table 3), but have led to a number 
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of related conclusions which are described below. 

One recent analytical study of the polariza- 
tion properties of masers by Elitzur (1998) has 
suggested that the usual thermal Zeeman formu- 
lation, V = ZC\ B \ dl/dv with C a constant, 
may not be valid for unsaturated masers. Unfor- 
tunately, a key uncertainty in the study of masers 
is that their degree of saturation cannot be ob- 
served directly. However, Elitzur (1998) also pro- 
poses that a maser's saturation state can be re- 
vealed by the ratio of its Stokes V profile to the 
derivative of Stokes I (i.e. K=Y /{dl/dv)). This 
is a consequence of the fact that maser amplifi- 
cation during unsaturated maser growth is expo- 
nential, causing a narrowing of the maser line. 
Under these circumstances, Elitzur (1998) pre- 
dicts that TZ= ZC\ B | is not a constant (i.e. 
C / constant) and instead assumes a Gaussian 
absorption shape. 

A test of this prediction was performed for the 
TZ profiles of W51C feature 2 (dynamic range = 
720:1) and the highest dynamic range maser in 
the Claussen et al. (1999) VLBA study of W28 
(W28 F39 [A]; dynamic range = 60:1). Note, 
that the VLBA OH maser data have three times 
the spectral resolution of the VLA data. The 
TZ profiles for both of these masers are relatively 
flat (to within their rms) and show no evidence 
of the Gaussian absorption shape predicted for 
unsaturated maser emission. We conclude how- 
ever, that while these TZ profiles are indicative 
of saturated maser emission, this test is some- 
what impractical for real data. The derivative 
of I passes through zero at line center by defini- 
tion, and toward the line wings as well depending 
on the rms of the spectra. A more useful varia- 
tion of this test is to investigate the linear depen- 
dence of Stokes V as a function of dl/dv since this 
test does not suffer from the discontinuities that 
are inevitable with the TZ test. The resulting V 
vs. dl/dv plots reveal the straight lines expected 
from saturated emission (i.e. ZC\ B |= con- 
stant) for W51C OH(2) and W28 F39 [A] along 
with several more of the masers in this study (i.e. 
G349.7+0.2 OH(3), CTB37A OH(l-5)), but re- 
main inconclusive for the lowest S/N masers (i.e. 
CTB33 OH(l)). 



Another perspective on a maser's degree of 
saturation can be gained by comparing an upper 
limit for its optical depth (T max ) assuming that 
it is unsaturated, with a lower limit (r m j n ) based 
on the flatness of its observed TZ ratio. For unsat- 
urated masers, T b = T x exp(T max ) (where T max is 
the maser's optical depth at line center and T x is 
the excitation temperature; see Elitzur 1998 for 
details). As mentioned above, the lower limit de- 
pends on the confidence with which a flat TZ pro- 
file can be determined, so that T m i n — e 
where e is the dynamic range of the observa- 
tions and V is a Bg /B e (Elitzur 1998). Elitzur 
(1998) notes that whenever these two limits are 
inconsistent (r m j n > T max ), the maser must be 
saturated. This 'inconsistency' was found to be 
true for every maser in our sample, indicating 
that OH (1720 MHz) masers may indeed be sat- 
urated. This was also found to be the case in 
a OH (1720 MHz) maser study of the galactic 
center by Yusef-Zadeh et al. (1999a). We point 
out, however, that we have only been able to cal- 
culate a lower limit for T b since the maser spots 
are not resolved at the present resolutions. This 
means that if the true values of (and hence 
T max) are large enough, the apparent saturated 
maser inconsistencies (r m j ra > T max ) could dis- 
appear. However, the sizes found for OH (1720 
MHz) masers in the W28 and W44 VLBA obser- 
vations of Claussen et al. (1999) ~ 50 — 180 mas, 
if typical, imply that the majority of masers in 
this study would still have T m i n > T max as re- 
quired for saturated maser emission. 

Even if these masers are likely to be satu- 
rated, a second intrinsic uncertainty in perform- 
ing Zeeman analysis on OH masers is the value 
that should be taken for the constant "C" in 
V= ZC\ B 1 1'. According to Elitzur (1998), 
this constant is modified from its thermal value: 
C t h = cos 6 to C m = 8/{3pcos 6) for saturated 
masers, where 6 is the angle between the mag- 
netic field vector and the line-of-sight. This so- 
lution to the maser polarization problem takes 
into account the fact that a photon generated 
via stimulated emission does not necessarily have 
the same polarization as the parent photon, and 
that the radiation is beamed. The parameter p in 
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this formulation depends on the geometry of the 
masing region, with p = 1 or 2 for filamentary or 
planar geometry respectively. Notice that in ad- 
dition to the factor of 8/(3p) difference between 
the thermal and maser equations for | B |, they 
also have a completely different dependence on 
cos 6. That is, while thermal radiation samples 
Bi os =\ B | cos 6, masers sample -B m oc| B \ /cos 6 
according to the Elitzur model. Therefore, this 
model suggests that the magnetic field needed to 
produce a maser's Stokes V signal is less than the 
field needed to produce a comparable signal from 
thermal radiation (i.e. C m > C t h)- For example, 
using | B | = (B e /C rn ), the Elitzur (1998) model 
predicts | B \ = 0.3pB e for 9 = 35°, and | B \ = 
O.lpBo for 6 = 75°. 

Ideally, the Stokes parameters U, Q, and V 
should only depend on 9 and | B |. Therefore, 
if analytic expressions for the Stokes parameters 
are known, it should be possible to calculate 6 
and | B | explicitly. According to Elitzur (1996), 
the analytic expression for the linear polarization 
fraction (| q | = (Q 2 + U 2 ) 1 / 2 /I) has the approxi- 
mate form 

2 



1 - 



3 sin 2 9 



(1) 



as long as xb « 1- Therefore, if the linear po- 
larization fraction is known, Eq. [1] can provide 
an estimate for 9 and hence the magnitude of 
C m . For the case of W51C OH(2), where we have 
measured q = 3.5% (§3.2.1), we obtain 9 rs 56° 
and C m ~ 5 /p. Consequently, in the Elitzur 
model, W51C OH(2) has | B |« (0.4 x p) mG 
compared to its estimated thermal Bq value of 
1.9 mG (where p = 1 or 2). 

Unfortunately, this technique can only pro- 
duce accurate results if the "total" linear polar- 
ized intensity is measured. The more extensive 
Stokes Q and U observations of OH (1720 MHz) 
masers in W28 and W44 (Claussen et al. 1997), 
indicate that q only ranges from 0.5% to 22% for 
all 30 of the masers with positive Q and U de- 
tections. In fact, the linear polarization of the 
masers in W28 and W44 is quite low with aver- 
age values of 4 ± 5% and 10 ± 5% respectively , 
similar to that found for W51C OH(2) (± gives 
the standard deviation of the distribution). It 



would be quite remarkable if q, and hence the 
angle between the line-of-sight and the magnetic 
field (9), were so similar in three different SNRs. 
Indeed it would be more natural for the polar- 
ized intensity of these masers to be more widely 
varied as is the case for SiO masers (see for ex- 
ample Mcintosh & Predmore 1993). There are at 
least two other effects in addition to the viewing 
angle (9) which can decrease the magnitude of q. 
These are 1) Faraday depolarization within the 
masing region and 2) curvature of the magnetic 
field lines within the masing region causing can- 
cellation within the beam (i.e. MHD waves, see 
Elitzur 1992). 

Faraday rotation can only be a significant fac- 
tor if the Faraday rotation sizescale is smaller 
than the length scale of the masing region (see 
Elitzur 1992 for details). From Elitzur (1992) the 
Faraday rotation sizescale is lp = 2xl0 16 (A 2 n e i3; os )~ 1 
where A is the observing frequency in cm, n e 
is the electron density in cm" 3 , and Bi os is the 
magnitude of the magnetic field along the line 
of sight in Gauss. Theoretical modeling of OH 
(1720 MHz) masers by Lockett et al. (1999) and 
Wardle (1999) suggest that the ionization frac- 
tion needed to produced a large enough column 
of OH for strong maser action is in the range 
10 -7 & Xe ^ 10~ 5 . These authors also find that 
the density in the masing region must lie in the 
range 1 x 10 4 cm -3 ^ <J 5 x 10 5 cm" 3 . Given 
our average magnetic field detection of Bq~ 1 
mG, the Faraday rotation length scale can be 
written 



l F = 6.2x10 
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/ n H 2 \ ( Xe \ 

UO 5 cm- 3 / UO- 6 / 



B, 



las 



1 mG 



cm. 



(2) 

This lp estimate is only three times larger than 
the OH (1720 MHz) maser gain length estimated 
by Lockett et al. (1999) of ~ 2 x 10 17 cm, imply- 
ing that Faraday depolarization may contribute 
to the low linear polarization intensities observed 
by Claussen et al (1997). The decrease of polar- 
ization intensity with increasing J in SiO masers 
is also thought to be a consequence of such Fara- 
day depolarization (Mcintosh & Predmore 1993; 
or Wallin & Watson 1997 for an opposing view). 
Unfortunately, the magnitude of such depolar- 
ization is difficult to quantify and quite model 
dependent (see e.g. Wallin & Watson 1997). 
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In addition, some evidence of tangling in the 
magnetic field lines may be indicated from the 
ten-fold increase in Bg measured toward W28 
F39 [A] with the VLBA (~ 2 mG) compared 
to the average value measured throughout W28 
with the VLA (~ 0.2 mG; Claussen et al. 1999). 
However, this particular spot did not have a pos- 
itive VLA Bg detection (due to the presence of 
multiple spatial components at the lower resolu- 
tion) so a direct comparison is not possible. More 
VLA vs. VLBA magnetic field strength compar- 
isons are needed to determine whether this ap- 
parent increase in Bg with higher resolution is 
real and/or common. 

These examples suggest that the linear polar- 
ization intensities observed toward W51C, W28, 
and W44 could be reduced from their 'true' val- 
ues by one or both of these effects, making the 
calculation of 6 directly from q unreliable. Al- 
ternatively, using the condition that q 2 + v 2 < 1, 
(where v =V/I), it is possible to obtain the fol- 
lowing constraint on 8 at the frequency of the 
peak in Stokes V, 16(xb) 2 <J cos 2 6 <; 2/3 as 
long as xb < 0.2 (see Elitzur 1996 for details). 
This means that the angle between B and the 
line of sight must be greater than ~ 35° in or- 
der for polarized emission to be observed at all 
(see also Goldreich, Keeley, &: Kwan 1973). This 
lower limit on (~ 35°), suggests that the mag- 
netic field values reported in Table 3 (Bg) are 
overestimated by factors of at least 1.5 (p = 2). 
The upper bound on 9 cannot be accurately cal- 
culated from these data since it depends on xb 
[ZBgj Av£>) and our current spectral resolution 
(0.27 km s _1 ; the highest available with the VLA 
in 2IF mode) is insufficient to resolve line split- 
tings with xb < 0.2. However, since Bg is an up- 
per limit on the total field strength (| B |), a use- 
ful limit on xb can be made by taking the average 
value of Bg/ Avfwhm using data from Tables 2 
and 3, in which case xb <^ 0.12 on the average 
and 6 <; 61°. This upper limit to 9 suggests that 
the average total field strength is | B |^ 0.2pBg, 
i.e. the thermal approximation (Bg in Table 3) 
overestimates the maser field ( | B | ) by less than 
a factor of C m ~ 5 /p. 

It is important to note, however, that the 



analytic expressions formulated by Elitzur are 
controversial and that other maser studies have 
reached different conclusions. For example, the 
numerical simulations of Nedoluha & Watson 
(1992) suggest that if xb << 1 and the maser is 
not strongly saturated, the thermal Zeeman rela- 
tionship is a good approximation (i.e. C m ~ 1). 
Unfortunately, it is quite difficult to compare the 
two methods since the Elitzur model is an an- 
alytical treatment, while the Nedoluha &; Wat- 
son (1992) result arises from a numerical simu- 
lation. Also notable is the Nedoluha & Watson 
(1990) finding that it is very difficult to produce 
linear polarization in masers unless they are at 
least partially saturated. It is impossible to make 
any quantitative statements on this issue with 
the data presented here since we only attempted 
linear polarization measurements toward one of 
our sources (W51C). However, the Claussen et 
al. (1997) OH (1720 MHz) maser study of all 
four Stokes parameters toward the SNRs W28, 
W44, and IC443 contains a total of 49 masers 
with S v > 200 mJy. About 60% of these have 
positive Q and/or U detections. In addition, all 
13 of the Claussen et al. (1997) masers with 3cr 
Bg detections also have positive Stokes Q and/or 
U detections. These comparisons suggest that in 
the Nedoluha & Watson model strong OH (1720 
MHz) masers are also likely to be at least par- 
tially saturated. 

An additional observational test offered by 
maser theory is to check whether there is any 
correlation between the linear polarization po- 
sition angle (P. A.) and the SNR shock front in 
the vicinity of a maser spot. Several studies have 
suggested that OH (1720 MHz) masers must arise 
in shocks which are propagating transverse to the 
line of sight (Lockett et al. 1999; Wardle 1999; 
Frail & Mitchell 1998; Claussen et al. 1997). It 
is also likely that the magnetic field vector lies 
preferentially in the plane parallel to the shock 
front since only this component can be amplified 
by shock compression. In addition, the linear po- 
larization P.A. for masers, unlike that of thermal 
dust emission, can be either parallel or perpen- 
dicular to the direction of the magnetic field in 
the plane of the sky. This difference occurs be- 
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cause the asymmetry of dust grains removes one 
degree of freedom. For these reasons, we might 
expect the observed P.A. of the linear polariza- 
tion to be either parallel or perpendicular to the 
shock front. Of course such a comparison is only 
meaningful if significant Faraday rotation of the 
P.A. has not occurred either because the Fara- 
day depolarization is insignificant or because the 
P.A. of the linear polarization is unaffected by 
it. A lack of net P.A. rotation in the presence of 
depolarization is suggested by Elitzur (1992) and 
the SiO maser observations of Mcintosh & Pred- 
more (1993), although Wallin & Watson (1997) 
present an alternate viewpoint based on numeri- 
cal simulations. 

For the case of W44 region E, the linear po- 
larization angle is ~ —28° (Claussen et al. 1997) 
and the OH (1720 MHz) masers are also dis- 
tributed along a line which is oriented NW/SE 
at about the same angle. Further indication of 
the orientation of the shock front comes from CO 
(J = 3 — 2) observations by Frail & Mitchell 
(1998) which show a ridge of shocked CO emis- 
sion parallel to the line of masers. This example 
suggests that there is a correlation between OH 
(1720 MHz) linear polarization position angles 
and the orientation of the shock front. The pres- 
ence of such a correlation could argue against the 
possibility of significant Faraday rotation within 
the SNRs or may result from Faraday depolariza- 
tion without significant net rotation. In any case, 
additional linear polarization measurements in 
conjunction with molecular observations of OH 
(1720 MHz) maser regions are needed to confirm 
this connection. 

These comparisons between observation and 
theory of OH (1720 MHz) masers have led to the 
following conclusions: (1) It is likely that OH 
(1720 MHz) masers are saturated based on their 
flat 1Z profiles (determined directly and from es- 
timates of maser optical depths; Elitzur 1998). 
This conclusion is also supported by the high 
incidence of linear polarization in these masers 
(Nedoluha & Watson 1992; Claussen et al. 1997). 
(2) The linear polarization intensities of OH 
(1720 MHz) masers (needed to calculate | B | in 
the Elitzur model) may be significantly reduced 



from their intrinsic values by Faraday rotation 
and/or by tangling in the magnetic field lines. 
This also means that it is impossible to calculate 
the reduction of | B | compared to Bq using these 
data directly (| B | oc BqcosO). However, for 
angles between 35° and 75°, the Elitzur model 
predicts that the total magnetic field strength 
should range from | B \ = (0.3 — 0.1)pBg (where 
p = 1 or 2 for filamentary or planar geometry re- 
spectively). Alternatively the numerical simula- 
tions of Nedoluha & Watson (1992) suggest that 
this correction factor is essentially unity. (3) It 
may be possible to predict the direction of the 
shock in which an OH (1720 MHz) maser arises 
from the position angle of its observed linear po- 
larization (as long as the P.A. is not significantly 
rotated by Faraday rotation). 

4.2. Magnetic Fields in Shocked Molec- 
ular Gas 

In this study we have reported ten new mea- 
surements of OH (1720 MHz) Zeeman magnetic 
field strengths in five galactic SNRs. Previous 
studies of this type (Claussen et al. 1997; Ko- 
ralesky et al. 1998; Yusef-Zadeh et al. 1996) 
have measured Bq in an additional five SNRs. 
The magnitude of Bq in all of these SNRs (in- 
cluding those measured here) have typical values 
between 0.1 and 4 mG. However, as the discus- 
sion in §4.1 demonstrates, the conversion factor 
between thermal estimates of Bq and the true 
maser field | B | is uncertain. Therefore, al- 
though we will continue to reference our mea- 
sured fields Bq (Table 3), it should be kept in 
mind that Bq could be an upper limit that over- 
estimates the true field by less than a factor of 
five. 

Maser theory suggests that OH (1720 MHz) 
masers originate in shocked molecular clouds 
(e.g. Lockett et al. 1999; Wardle 1999). If 
so, the measurements reported here of Bq (and 
in the references cited above) must originate in 
post-shock gas. These theories suggest that OH 
(1720 MHz) masers can only be pumped effi- 
ciently for densities in the range 1 x 10 4 cm -3 
<J uh 2 <; 5 x 10 5 cm~ 3 and temperatures in the 
range 50 K ^ T ^ 125 K (see Lockett et al. 
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1999). Indeed, when independent measurements 
of the conditions in the post-shock gas have been 
made (see Frail & Mitchell; Reach k Rho 1998, 
1999) the gas properties are in agreement with 
these theoretical expectations. 

It remains an open question how magnetic 
fields of these strengths are generated; i.e. shock 
compression vs. turbulent amplification (Jun & 
Norman 1996). In what follows we will show 
that compression of the existing ambient molec- 
ular cloud field is all that is required to produce 
the observed field strengths. One further argu- 
ment against significant turbulent amplification 
of the fields is the likelihood of destroying the 
maser action due to loss of velocity coherence in 
a turbulent velocity field. 

A number of Zeeman studies have been under- 
taken in the past decade to detect magnetic fields 
in molecular clouds for the purpose of studying 
star formation (e.g. Roberts et al. 1993; Bro- 
gan et al. 1999). In a recent review by Crutcher 
(1999) of the existing data for star forming re- 
gions, the magnetic field was found to scale with 
density as | B |oc n 0A7 . Crutcher notes that 
there are two possible physical interpretations 
for this relationship: (1) Such a relationship be- 
tween | B | and n has been predicted by Fiedler & 
Mouschovias (1993) based on studies of ambipo- 
lar diffusion; (2) A similar relation is suggested 
by the observed invariance of the Alfvenic Mach 
number vtia = V3o~/Va ~ 1 in molecular clouds, 
where V A =\ B \ /Airp 1 / 2 and a = Av/(8 In 2) 1 / 2 
(see Bertoldi & McKee 1992; Zweibel & McKee 
1995). This invariance implies | B |oc Av^fp. 

For the range of gas densities expected in 
OH (1720 MHz) maser regions (1 x 10 4 cm~ 3 
£ n H2 <, 5 x 10 5 cm" 3 ; Lockett et al. 1999), 
the range of magnetic fields predicted by the re- 
lation | B |oc n - 47 is 75 pG to 475 pG. Clearly 
the OH (1720 MHz) Bq measurements reported 
in this work (Table 3) are greater than those 
predicted by Crutcher's relation (see also Figure 
14). However, there is a great deal of uncertainty 
in this statement since Crutcher considers B\ os 
(Bios=\ B | cos 6) which underestimates the field 
by factors of ~ 2 while our determination of Bg 
could overestimate the field by as much as a fac- 
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Fig. 14. — Plot of Bgot n aA7 (Crutcher 1999; solid 
thick line) with the OH (1720 MHz) mascr's Bg pa- 
rameter space superposed. The two dashed lines show 
the lcr errors on the fit obtained by Crutcher. 

tor of five according to the Elitzur model (see 
§4.1). 

If we assume, however, that the two types 
of magnetic field measurements (molecular cloud 
vs. OH (1720 MHz) masers) are approximately 
comparable (i.e. the thermal Zeeman equation 
is valid for masers) it may not be surprising 
that their magnitudes are different. This is be- 
cause theory and observations suggest that OH 
(1720 MHz) masers do not originate in undis- 
turbed "normal" molecular clouds, but rather 
have experienced a shock (see above). If this 
is the case, the Bertoldi &: McKee interpreta- 
tion for the scaling of | B | (oc Av^fp) suggests 
that OH (1720 MHz) masing regions simply have 
larger linewidths for a given density than un- 
shocked molecular clouds. If we set Va = V3o 
as suggested by Bertoldi and McKee (1992), then 
| B |= OAAvrip where n p is the proton density 
in cm -3 and the line width (Av) is in km s _1 . For 
the case of the W51C OH (1720 MHz) masers, 
we can take Av = 10 km s~ 4 based on the CO 
observations of Koo (1999), and n p = 1 x 10 5 
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cm 3 from the typical density needed to excite 
OH (1720 MHz) masers (Lockett et al. 1999; 
Wardle et al. 1999). With these estimates, the 
implied magnetic field strength is | B | = 1.3 mG, 
in close agreement with the Bg values observed 
toward W51C. 

Alternatively, if the Bg's measured in OH 
(1720 MHz) masers should be reduced by factors 
of <; 5 (i.e. Elitzur's maser polarization model), 
the fields estimated from OH (1720 MHz) masers 
are in agreement with the values predicted by 
Crutcher's relation. In this case, the magnetic 
fields in masing regions and molecular clouds fol- 
low the same scaling with density without shock 
amplification of the field, i.e. no Av dependence. 
Such agreement might be expected if the ob- 
served scaling were the result of ambipolar dif- 
fusion in both types of regions. Clearly con- 
vergence on our understanding of the nature of 
maser polarization (and hence | B |) is needed to 
distinguish between the two possibilities. 

In a radiative shock, compression of the gas 
follows the relation r] = \^2V s /Va , where Va 
is the Alfven velocity in the pre-shock gas (see 
Draine & McKee 1993). Using the relation 

1/2 

B oc n , it can be shown that Va — 2 
km s _1 (see also Heiles et al. 1993). Thus 
for reasonable shock velocities (10 - 50 km s -1 ; 
Lockett et al. 1999; Frail & Mitchell 1998), 
rf = 7—35. This compression is sufficient to 
enhance B Q to Bg without invoking turbulent en- 
hancement of the magnetic field, assuming that 
B ps /B = 7] (see Chevalier 1999; Frail & Mitchell 
1998). For example, if we take a typical shock 
speed of V s ~ 25 km s _1 and typical values 
of B Q ~ 70 /j,G for molecular clouds of density 
n Q ~ 5 x 10 3 cm" 3 (Frail & Mitchell 1998; Wardle 
1999; Crutcher 1999), 77 » 18 and B ps ~ 1.2 mG 
- in close agreement with the observed values of 
Bg. This demonstration shows that typical cloud 
and shock parameters can lead to enhancement 
of the pre-shock magnetic fields to the observed 
values of Bg. 



4.3. Energetics Implied by Observed 

Bg 

Using the thermal Zeeman equation, we have 
detected magnetic fields (Bg) between 0.2 and 
2 mG (Table 3) in OH (1720 MHz) masers as- 
sociated with five galactic SNRs. However, as 
described in §4.1 these thermal field strengths 
may overestimate the "true" field strength | B \ 
by a factor between ~ 2 — 5 (i.e Elitzur 1998). 
Alternatively, other theoretical studies suggest 
that Bg is a good approximation to | B | (i.e. 
Nedoluha & Watson 1992). For this reason, along 
with our inability to calculate the exact correc- 
tion factor in the Elitzur model, we have chosen 
to use Bg in the calculations below. However, 
the fact that Bg could be an upper limit to the 
magnetic field strength should be kept in mind. 

Regardless of the exact details of the Bg to 
I B I conversion or the origin of these fields (ar- 
gued to be shock compression of the ambient field 
above), it is clear that the OH (1720 MHz) fields 
are strong (compared to those found in dark 
clouds ~ 30 /iG, for example). This is an impor- 
tant finding since the magnetic field plays a key 
role in many aspects of the SNR/molecular cloud 
interaction including : (1) As noted above, the 
magnetic field determines and, in fact, limits the 
amount of shock compression in the post-shock 
gas. (2) Magnetic support may help stabilize the 
post-shock cloud, via the magnetic force acting 
on ions perpendicular to the magnetic field lines. 
(3) In the model of Lockett et al. (1999), heating 
by ambipolar diffusion is needed to extend the 
length of time the post-shock gas spends at tem- 
peratures favorable to OH (1720 MHz) maser in- 
version. Ambipolar diffusion is also thought to be 
the method by which subcritical clouds dissipate 
their magnetic energy and form stars (Ciolek & 
Mouschovias 1995). Given the many roles that 
the magnetic field can play in the evolution of a 
shocked cloud, we estimate the magnetic energy 
compared to other energy sources within the re- 
gion below. 

From the values of Bg listed in Table 3 the 
magnetic pressure (Pb ps = Bp S /8ir), in the post- 
shock gas of the five SNRs studied here, ranges 
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from 1CT 7 to 1CT 9 erg cm . These values are 
large compared to both Pism (~ 5 x 10~ 13 erg 
cm~ 3 ) and Pthermai of the hot X-ray emitting gas 
interior to the remnant, ~ 2 x 10~ 10 erg cm -3 
(Kulkarni & Heiles 1988; Frail & Mitchell 1998; 
Claussen et al. 1997). From the X-ray observa- 
tions of Koo et al. (1995) of W51C (n e = 0.3 
cm -3 ; T e = 3 x 10 6 K), we can estimate directly 
that Pthermai = 2n e kT e « 3 x 1CT 10 erg cm -3 
for this SNR. The ram pressure of the shock can 
also be compared to Pg since they should be ap- 
proximately equal if the field strengths are pro- 
portional to the amount of shock compression, 
i.e. Bp S /8ir = p Vg. For example, using the val- 
ues of n and V s estimated in §4.2 (n = 5 x 10 3 
cur 3 and V s = 25 km s" 1 ), p Q V^ = 5 x 1(T 8 
erg cm -3 . The equivalent magnetic pressure for 
W51C OH(l) (B g = 1.5 mG) is 9x 10" 8 erg cm" 3 , 
in close agreement with the estimated ram pres- 
sure. Therefore, the magnetic pressure in the 
post-shock gas dominates over the thermal pres- 
sure of the SNR but is approximately equal to 
its ram pressure. The fact that the ram pressure 
is almost three orders of magnitude larger than 
the thermal pressure of the shock indicates that 
this SNR is in the radiative phase as expected for 
SNRs with OH (1720 MHz) masers. 

5. SUMMARY AND CONCLUSIONS 

We have observed the OH (1720 MHz) line in 
five galactic SNRs to measure their magnetic field 
strengths using the Zeeman effect. We detected 
all 12 of the bright (S u > 200 mJy) OH (1720 
MHz) masers previously observed by Frail et al. 
(1996) and Green et al. (1997) and measured 
significant magnetic fields (i.e. > 3a) in ten of 
them. The estimated field strengths, which we 
denote Bg, range from 0.2 to 2 mG and are in 
good agreement with those measured in the five 
other SNR (0.1-4 mG) for which Zeeman OH 
(1720 MHz) maser studies exist (Koralesky et al. 
1997 and references therein). In these studies, 
the field strengths were calculated using the ther- 
mal Zeeman equation for Zeeman splitting less 
<< the line width (i.e. V = ZBgdI/du). In this 
formula, the total field strength (| B |) is related 
to Bq by | B \=Bg/C m where C m is a function 



of the viewing angle 6 (C rn = cos 6 for thermal 
emission). However, there exists some contro- 
versy about the 6 dependence of C m for saturated 
masers and whether it is even a constant for un- 
saturated masers (see §4.1). For example, in the 
Elitzur (1998) model C m oc \/{cos6) for masers. 
From comparison of these data with maser polar- 
ization studies we conclude that these OH (1720 
MHz) masers are likely to be saturated and that 
C m ^ 5, i.e. the Bg values reported in Table 3, 
overestimate the actual field strengths (\ B \) by 
less than a factor of five. 

We estimate that these magnetic field strengths 
are consistent with the hypothesis that ambient 
molecular cloud B-fields are compressed via the 
SNR shock to the observed values (§4.2). Indeed, 
field strengths of this magnitude exert a consid- 
erable influence on the properties of the shocked 
cloud. For example the magnetic pressures es- 
timated from the values of Bg listed in Table 3 
(io- 7 -io~ 9 erg cm ) far exceed the pressure in 
the ISM or even the thermal pressure of the hot 
gas interior to the remnant (§4.3). We also find 
that this range of magnetic pressures is in very 
good agreement with the ram pressure expected 
from C-type radiative shocks. 

In §4.2 and §4.3 we show that there is excellent 
agreement between our values of Bg and those 
implied by shock compression and ram pressure 
using typical values from the literature. It is 
somewhat difficult to maintain this agreement if 
Bg should be reduced by as much SIS 3i factor 
of ^ 5 (i.e. Elitzur's maser polarization model 
§4.1), although there is probably sufficient un- 
certainty in the parameters to allow such fine- 
tuning. We also show that the observed values 
of Bg for OH (1720 MHz) masers are greater than 
those observed in molecular clouds for the same 
range of densities (Crutcher 1999). It is possible 
that this difference is due the intrinsic physical 
nature of the | B |oc n 1 / 2 relation (§4.2) or may 
be the result of overestimating the maser B-fields 
(§4-1). 

In the future, knowledge of the field strength 
could be used in conjunction with molecular data 
to study the physics of these molecular shocks in 
more detail. In one such study Frail & Mitchell 
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(1998) mapped the distribution of molecular gas 
in the vicinity of several masers spots in W28 
and W44 with the JCMT. These observations re- 
vealed that the OH masers are preferentially lo- 
cated along the edges of thin filaments or clumps 
of molecular gas, suggesting compression of the 
gas by the passing shock. In addition to this 
morphological evidence, the density, temperature 
and velocity dispersion of the gas at these loca- 
tions suggested that the OH (1720 MHz) masers 
originate in post-shock gas. Combining the VLA 
and JCMT data they were able to show directly 
that the magnetic pressure dominates over the 
thermal pressure in the post-shock gas, balanc- 
ing against the ram pressure of the gas enter- 
ing the shock wave (i.e. Bp S /87r = p V|). Thus 
OH (1720 MHz) measurements of Bg lead to con- 
straints on the physics of these molecular shocks 
that are difficult to obtain any other way. When 
used in conjunction with molecular observations, 
it should be possible to fully specify the proper- 
ties (i.e. geometry, density, temperature, veloc- 
ity) of C-type shocks. 
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